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(57) ABSTRACT 

A fast-converging, computationally simple, method for rec- 
ognizing a single frequency tone or a sinusoid in a signal 
without prior knowledge of the tone frequency. The method 
employs a second order or higher auto -regressive model and 
includes: (a) sampling the signal at a constant sampling rate, 
and, for each sample, recursively determining a finite num- 
ber of correlation coefficients using a time -reversed, expo- 
nentially weighted, future sliding equivalent of the signal, 
wherein the correlation coefficients are determined using 
pre-existing values of the correlation coefficients determined 
in a previous iteration, a current sample of the signal and at 
least two consecutively previous samples of the signal; (b) 
periodically determining at least the second auto -regressive 
coefficient modeling the signal using the correlation coeffi- 
cients; and (c) recognizing the presence of the tone based on 
the value of the second auto-regressive coefficient. 

18 Claims, 3 Drawing Sheets 
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RECOGNITION OF A SINGLE FREQUENCY 
TONE 

FIELD OF THE INVENTION 

The invention relates to the art of recognizing a single 
frequency tone in a signal, and more particularly, the rec- 
ognition of such a tone in a signal without any prior 
knowledge or expectation of the tone frequency. 

BACKGROUND OF THE INVENTION 

In echo-cancelers, it is necessary to freeze or disable the 
echo-cancelling operation when specific tones are propa- 
gated through the echo-canceller. One well-known case of 
such a tone is the AC15 tone used for inter-PBX signalling 
(as set forth in Signalling System AC15, British Telecom 
Network Requirement Document No. BTNR181). PBXs 
using this signalling scheme do so by using a single fre- 
quency tone of 2280 Hz. It is possible for two interconnected 
PBXs to generate the 2280 Hz tone simultaneously. Thus, if 
an echo -canceller is in the exchange path of the two PBXs, 
the echo-canceller could eliminate the tone in one direction 
simply because the tone in this direction can be considered 
as the echo of the signal in the opposite direction. 

There are many other standardized bi-directional signal 
tone exchanges. For example: 

1) Continuity testing (1004 Hz); 

2) AC15 (2280Hz); 

3) ITU-T Signalling System No. 4 (2040 Hz and 2400 
Hz); 

4) ITU-T Signalling System No. 5 (2400 Hz and 2600 
Hz); 

5) ITU-T Signalling System No. 6 and No. 7 (2010 Hz); 

6) China PIT TS-02 (2600 Hz. 
These, as well as many other standardized single tone 

exchanges, employ numerous different single tone frequen- 
cies which must be detected for various purposes. 

One way to solve the tone detection problem is to tailor 
multiple specific-tone detectors to detect every one of the 40 
frequencies involved. These prior art tone detectors are 
capable of detecting a single tone based on an expectation of 
what that tone should be. The drawback with these types of 
detectors is that their usability is restricted to the pre- 
programmed frequencies, i.e., one detector is required for 
each frequency one is trying to detect. Thus, the processing 
requirement grows linearly with the number of distinct 
frequencies, which must be detected. It would be more 
advantageous to have a more generic tone detector, whicb 
has the capability to recognize if a signal is a pure sinusoidal 
tone and, if, so, to determine the frequency thereof. Such a 
detector would have the distinct advantage of not having to 
know all the single frequencies involved in the network. 
Also, because such a detector can dynamically determine the 
frequency, various other benefits can be gained. For 55 
instance, it is possible to detect the echo protector tone used 
in the prologue of full-duplex voiceband data such as 
defined in V.25 EPT, V.8, or V.8bis standards (see, 
respectively, 

ITU-T Recommendation V.25 (1996), Automatic answer- 
ing equipment and general procedures for automatic calling 
equipment on the general switched telephone network 
including procedures for disabling of echo control devices 
for both manually and automatically established calls; 

ITU-T Recommendation V.8 (1994), Procedure for start- 
ing sessions of data transmission over the general switched 
telephone network; and 
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ITU-T Recommendation V.8bis (1996), Procedure for the 
identification and selection of common modes of operation 
between data circuit-terminating equipment (DCEs) and 
between data terminal equipment (DTEs) over the general 
switched telephone network and on leased point-to-point 
telephone- type circuits). 

A number of prior art techniques have been proposed to 
model sinusoidal signals. These include Prony's method or 
Pisarenko's spectral line decomposition. Unfortunately, 
these two methods use fairly complex mathematical opera- 
tions: Prony's method requires a polynomial root extraction 
and Pisarenko's technique requires calculating Eigenvalues 
and Eigenvectors. (See, respectively, deProny, Baron Gas- 
pard Riche, "Essai experimental et analytique: sur les lois de 
la dilatabilitc de fluides elastique et sur celles de la force 
expansive de la vapeur de ralkool, a differentes 
temperatures", Journal de Vfccole Poly 'technique 1795, 
volume 1, cahier 22, 24-76; and 

Pisarenko, V. F., "The Retrieval of Harmonics from a 
Covariance Function", Geophysical Journal of the Astro- 
nomical Society, 1973, 33:347-366). 

However, one technique, which tends to be mathemati- 
cally simpler, is the auto-regressive modeling or linear 
prediction coding. For instance, U.S. Pat. No. 5,495,526 to 
Cesaro et. al. and U.S. Pat. No. 5,619,565 to Cesaro et. al. 
disclose single frequency tone recognition techniques based 
on the auto-regressive model. The drawback with the tech- 
niques outlined in these patents is the imprecision of the 
algorithms used to determine the auto -regressive coeffi- 
cients. Secondly, these techniques use an adaptive filter, 
which consequently results in a comparatively long time 
required to detect a single tone frequency, i.e., poor conver- 
gence time. 

SUMMARY OF INVENTION 
The invention seeks to overcome many of the disadvan- 
tages of the prior art and more particularly provide a single 
frequency tone detector able to detect a tone in a signal 
without any prior knowledge or expectation of the tone 
frequency. 

Broadly speaking, one aspect of the invention provides a 
method for recognizing a tone in a signal using at least a 
second order non-biased auto -regressive model defined by a 
first auto-regressive coefficient and a second auto-regressive 
coefficient. The method includes: 

(a) sampling the signal at a constant sampling rate and, for 
each sample thereof, recursively determining a finite 
number of current correlation coef&cients using an 
exponentially weighted, future sliding equivalent of the 
signal. The signal is time reversed such that prior to 
each determination of the current correlation coeffi- 
cients as aforesaid, a current received sample and 
consecutively received previous samples of the signal 
are time reversed such that the current received sample 
is defined as an oldest sample, an oldest of the con- 
secutively received previous samples is defined as a 
current sample and, where more than two consecutively 
received previous samples of the signal are used in the 
determination, all the consecutively received previous 
samples other than the oldest thereof are also defined to 
be time reversed as to their received sampling order. 
The current correlation coefficients are determined 
using pre-existing values of the correlation coefficients 
determined in a previous determination, the current 
sample and at least two consecutively received previ- 
ous samples of the signal; 

(b) periodically determining at least a second auto- 
regressive coefficient for modeling the signal using the 
correlation coefficients; and 
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(c) recognizing the presence of the tone based on the preferred embodiments thereof and the accompanying draw- 
second auto -regressive coefficient. ings which illustrate the preferred embodiment by way of 

Another aspect of the invention provides a method for example only. In the drawings: 

detecting a tone and its frequency in a signal using at least FIG ^ is a ^^0^ block diagram of a single fre- 

a second order non-biased auto-regressive model defined by 5 qucn cy tone detector in accordance with the preferred 

a first auto-regressive coefficient and a second auto- embodiment; 

regressive coefficient. The method includes: u * * * ^ u *u * 

& . . , . FIG. IB is a flowchart of a program executed by the tone 

(a) sampling the signal at a constant sampling rate and, for detcctQr shown k nG ^ and 

each sample therof, recursively computing a finite . .... 

number of current correlation coefficients using an 10 A 2 » a s y s ' em block diagram showing the tone 

exponentially weighted future sliding equivalent of the detector <* 11145 preferred embodiment used to control an 

signal. The signal is time reversed such that prior to echo-canceler circuit. 

each computation of the current correlation coefficients DETAILED DESCRIPTION OF PREFERRED 

as aforesaid, a current received sample and consecu- EMBODIMENTS 

tively received previous samples of the signal are time 15 

reversed such that the current received sample is FIG. 1A shows a functional block diagram of a single 

defined as an oldest sample, an oldest of the consecu- frequency tone detector 10 according to the preferred 

tively received previous samples is defined as a current embodiment. The detector 10 is preferably implemented by 

sample and, where more than two consecutively a digital signal processor (DSP) 12 or other type of processor 

received previous samples of the signals are used in the 20 or hardware commercially available from a variety of 

computation, all the consecutively received previous sources. The DSP 12 comprises an input line 14 for receiv- 

samples other than the oldest thereof are also defined to ing an analogue signal, which is analyzed to determine if the 

be reversed as to their received sampling order. The signal currently exhibits the characteristics of a single tone 

current correlation coefficients are computed using frequency. The DSP 12 samples the signal on input line 14 

correlation coefficients computed in a previous 25 at a fixed sampling rate f, and, through a built-in (or 

iteration, the current sample of the signal and at least external) analogue-to-digital converter, provides digital val- 

two consecutively received previous samples of the ues indicative of the level of the signal. These samples are 

signal; used by an auto-regressive modelling or linear predictor 

(b) periodically computing at least the first auto- coding program executing on the DSP 12, as discussed 
regressive coefficient and the second auto-regressive 30 below, to detect the presence of tone. If a tone is detected, 
coefficient for modelling the signal using the current the DSP 12 sets a Boolean output detect signal 16 and 
correlation coefficients; optionally provides the value of the detected frequency on 

(c) recognizing the presence of the tone based on the an output frequency signal 18. When no tone is detected, the 
second auto -regressive coefficient; and detect signal 16 is reset and the frequency signal 18 is 

(d) determining the frequency of the tone based on the 35 mvalld - 

first auto-regressive coefficient and the sampling fre- As mentioned, the DSP 12 employs auto-regressive mod- 

quency. eling. The basic principle behind this type of modeling is to 

In the preferred embodiment of the invention, at least the consider the analyzed signal as being the output of a 

first two auto-regressive coefficients modeling the signal are recursive filter excited by white noise. This modeling 

periodically computed using the correlation coefficients, and 40 approach directly translates in the frequency domain as 

the frequency of the tone is determined based on the value trying to match the transfer function of the filter to the 

of the first auto-regressive coefficient and the sampling spectral contents of the analyzed signal. The transfer func- 

frequency. uon °f a V th order auto -regressive model is: / 

Moreover, the step of periodically determining at least a 

second auto-regressive coefficient is executed every P 45 H{z)= - ^ 

samples, P being an integer number of at least one. l + aiZrl + W 2 + .»+V* 

Moreover, the step of recognizing the tone preferably 

includes confirming that the second auto-regressive coeffi- ^ ^ afe ^ aut ressive model coefficients and 

cient is within a predetermined range of foe square of an js a ^ ^ ^ auto . r cssive model has its ^ 

exponent.^ decay factor used to define the exponentially £ ^ ^ ^ ^ where ^ H(z) 

weighted futare shdmg equivalent of the Sl gnal. The pre- Action becomes a linear prediction, as follows: 
ferred embodiment additionally validates the presence of the 

tone by determining the power of the signal and confirming *„=pw„-ai*„-i-«2*„-2- • • • ( 2 ) 
that the power exceeds a threshold power, and, optionally, by 

computing the first auto-regressive coefficient and confirm- 5S where p is a gain factor and u„ is a prediction error of the 

ing the stability of the first auto-regressive coefficient over system (white noise). This equation indicates that the sample 

a number of iterations. ' x » can be predicted from the p previous samples with some 

It will be noted from the foregoing that no prior knowl- error: 
edge of the frequency of the detected tone is required. It 

should also be noted that a history of only three samples of 60 *„ = x n - pu„ « -«! - a 2 x„. 1 - ... - apx^ (3) 
the signal are required to implement the method and this, in 
conjunction with its recursive nature, enables the method to 

converge relatively quickly. Minimizing the prediction error un is achieved by resolv- 
ing the Yule-Walker set of equations. To detect a sinusoidal 

BRIEF DESCRIPTION OF DRAWINGS 65 f t ^ on i y necessary to use a second order model. This 

The foregoing and other aspects of the invention will can be understood by considering the representation of the 

become more apparent from the following description of the sinusoidal signal in the Z domain: the Z transform of a pure 
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sinusoid is represented by complex conjugate poles located debouncing to add robustness to the tone recognition. This 
exactly on the unit circle. debouncing operation will increase the time required for the 

The Yule Walker set of equations for a 2 nd order predictor detection, but minimizes the possibility of false detection. 
(p=2) is: A preferred way to overcome the time resolution limita- 

5 tion of the block processing approach is to perform the 
(4) computation of the a a and 04 coefficients at every new 
incoming sample. This in turn requires new correlation 
coefficients to be obtained at every sample. If the block 
processing approach is still retained, this technique has often 
10 been termed a "sliding window''. Using a sliding window 
where the a,- are the auto-regressive coefficients and the a,-,- involves executing a large number of computations at every 
are the correlation coefficients of the signal x„. If the gain incoming sample. This is undesired. In the preferred 
factor is not required or used, the previous set of equations embodiment a sliding exponential window is employed 
can be reduced to: instead of a rectangular window in order to minimize the 

15 computational requirements. Modeling an exponentially 

Ia n anirai] \ a n] (5) windowed sinusoid will not deliver the same o^ and a 2 

*2i *22 JU2 J = L«2o J coefficients for the auto-regressive modeling as would a 

normal sinusoid, but the effects will be explained in greater 
detail below. 

The characteristics of interest of a sinusoidal signal are 2Q Applying an exponential window on a signal means 
completely defined from the a 2 and a^ coefficients, as multiplying the samples by P", where p is a decay factor in 
discussed in greater detail below. ^ range 0 f 0 to 1. As n grows (toward the most recent 

One estimate R, y of the correlation coefficients a f/ of the samp i e ), the more the corresponding sample is attenuated. In 
signal can be obtained from the biased correlation esumation order t0 minimize computational requirements for the 
defined by: 25 correlation, a recursive computation is desired, i.e., where 

the correlation coefficients associated with the newest 
1 ^ ^ sample are computed from correlation coefficients associ- 

R v = jj/_^ n -M<n- ft atec j me p rev i ous sample. Therefore, the correlation 

computation according to the preferred embodiment 
30 employs a reverse exponential window where the most 
where x(n) is the n th sample of the signal and N is the recent sample is not attenuated, but the older the sample, the 
number of samples used to estimate the correlation coeffi- more attenuated it is. This is preferably accomplished by a 
cients. This estimate of the correlation implies that the data time-reversing operation. There is an undesired side effect in 
outside the N sample analysis window have a value of 0. attempting to time-reverse an exponentially decaying sinu- 
This creates a side effect where the auto-regressive model of 35 S oid: the signal to model then represents an unstable system 
a time-limited sinusoid is dependent on the initial and final and, technically speaking, auto-regressive modeling is not 
phase of the sinusoid, as well as the number of samples used valid on unstable signals. Time-reversing is not an issue 
to estimate the correlation coefficients (i.e., the "rectangular when using the biased correlation estimate [equation (6)] 
window" size). This can be better understood by noticing because the R 0 are the same as but this is not true for 

that the sinusoid is forced to be time limited and the Z 40 the other form [equation (7)] of the correlation estimate 
transform of this signal is no longer represented by a pair of where the a u are not the same as a^ +Jt . To avoid the 
complex conjugate poles located exactly on the unit circle. undesired side effect of the time reversing, the preferred 
A more preferred alternative to the biased correlation embodiment considers the incoming samples to be time 
estimation provided by equation (6) is to artificially reduce reversed. In other words, the most recent sample is always 
the number of samples used in the cross product [x(n-i)*x 45 tagged with the index 0. When a new sample needs to be 
(n-j)] by making sure that no sample is assumed to be zero. analyzed, the previous sample, originally at index 0, is 
This correlation function is denned by: considered being at index 1 and the new sample inherits the 

index 0. This means that the current sample is considered to 
j w-i (7) be the oldest sample, and the oldest sample is considered to 

^J- jj^Yi X ^ n ^^ n ~^ 50 be the current sample for the purpose of stabilizing the 

"~ p time-reversed signal. 

To understand how the correlation coefficient estimates 
where x(n) is the n'* sample of the signal, N is the number <*/, can be recursively computed according to the preferred 
of samples used to estimate the correlation coefficients and embodiment using a time-reversed exponentially weighted 
p is the order of the auto -regressive model. This estimate of 55 future sliding signal, consider an expansion of the compu- 
the correlation function does not suffer from the problems tation for all six a iy 
associated with the biased correlation estimate provided by 

equation (6). (omitting the scaling factor — — 1 

In the operation of any detector it is generally desirable to V n - pj 

detect the presence of a tone and optionally to determine its 60 

frequency as fast as possible. Conventionally, a block pro- used m ^ re(hlced second . order Yule-Walker set of equa- 

cessing approach has been used in performing the auto- ^ for &e Qrigmal signa] spanning the mdex rangc of ^xo 

regressive model; i.e., a number N of signal samples are tQ mfinit ^ zerQ ^ the mogt recent of ^ time 

accumulated and then processed. This approach implicitly reV ersed sequence: 

limits the detection resolution to twice the window size 65 

because the first analysis window may not constitute 'all a l0 -x 1 x 2 +^3+x 3 jf 4 + 

samples of a sinusoid. Moreover, there is a need for some a 20 =A^r 2 +x 1 x 3 +x 2 x 4 + 
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a il =x i x l +x^c 2 +x^+ where: 



a 1 ^x 1 xb+* 2 x 1 +x 9 x 2 + nf> (2nf 0 \ (13) 

a 22 -jc 0 x 0 +x 1 x 1 +x 2 x 2 + (8) 



ai= -2/fcos(^) and 
a^/J 2 (14) 



With the newest sample being indexed at -1 (index -1 is and 
used here to show the effect of the newest sample) this set 

of equation becomes: 10 Knowing the decay factor used for the correlation 

computation, it is possible to know when the analyzed signal 
a'io°^i+*i*2+*2*3+*>*4+ ^ a P ure SUTU soidal t° ne by verifying the computed value of 

a^ which equals to the square of the decay factor. If the 
a 2<r*-i*i +*b*2+*i*3+*2*4+ frequency of the tone is desired, the coefficient a 2 provides 

15 a value representing or indicative of the frequency, and if the 
a , 11 ^ o x 0 +x 1 x 1 +x 2 x 2 +^jf 3 + exact numerical value of the frequency f 0 must be reported 

it can be derived from: 



a.' 2X <*x_ x x 0 +xjc x +x 1 x 2 +x&+ 20 
az2~*-i*^ + *o*o + *i*i + *2*2 + ( 9 ) 



U „ l( « M (15) 



The first two a a and 2^ auto-regressive coefficients are 

It can be noticed that only the first element on the right side obtained by solving the second-order Yule-Walker set of 

of the set of equations (9) contributes to update the value of equations. If it is desired to merely recognize if a signal 

the new correlation coefficients. Dealing with an exponen- under analysis is sinusoidal, the value of the a 2 coefficient in 

tially attenuated sliding window, the sample x„ can be and of itself indicates the sinusoidal characteristics of a 

re-defined as P"x„ with the new x„ being the un-windowed signal. The coefficient ^ is obtained from the correlation 

original signal. The correlation coefficients can then be coefficients by the following equation: 
recursively updated by the following equations, considering 
x 0 to be the most recent sample, x a the previous sample, and 
x 2 the sample previous to Xj: 



30 



ff2lQ'L0-O'iiQr2o (16) 

a 2 = • 



<*'io-P 3 *i*2 + 0 2a io The coefficient a 17 which is indicative of the tone 

35 frequency, is obtained from the correlation coefficients and 

a'2o=P 2 *o*2 + P 2a 2o the a 2 coefficient as follows: 

a'n=P 2 *i*i+P 2 au aiQ + ai2a2 (17A) 

Qi = . 

40 

a' 21 -^x 0 x 1 +^ 2 a 21 Alternatively, the following equation may be used in view of 

equation (14): 

a22-*o*o+P 2 a22 ( 10 ) 

Thus, to recursively compute the correlation coefficients, a 45 a * - 

history of only the three most recent samples (x 0 , Xj and x^) 
is required as well as the previous values of the correlation 
coefficients. When updating the correlation coefficients, the preferred 

As mentioned above, the characteristics of a sinusoidal embodiment enables some reduction in the computation 

signal are completely defined from the a, and a* auto- 50 complexity to be achieved. First, it will be noticed that a, 2 
regressive coefficients. To understand why, consider an * ""f a * '°2i therefore only one update computation is 

0 . • 11 1 • . A n , .... required for these two correlation coemcients. Second, the 

exponent,ally decaying sme wave defined by: ^ Qf ^ M use& ^ fe p2 ^ ^ 

value of 0*22 which is using x 0 . In the same line, the current 
x n = Ay^sinU) + — n) 5S value for a 10 which uses x 1 and x^ is p 2 times the previous 

\ f* J value a 12 which is using Xq and x r Accordingly, to minimize 

the computations, the following set of equations can be 
performed in this order to recursively obtain the correlation 
where A is the sine wave amplitude, p is the decay factor, <P 0 coefficients: 
is the phase at time zero, f 0 is the frequency of the sine wave, 
and f, is the sampling frequency. This sine wave can be 60 <x'io-P 2a i2 
described using the auto-regressive model with: a . ^ 



A (12) 



a! 2Q =$ 1 x& 2 +& l au 1Q 



a'22^o^) + P 2 a22 
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From the foregoing, it will be noted that it is possible to and alternatively, 
blindly recognize a pure tone by using the a. 1 and a^ 

auto-regressive coefficients of a second order auto- \— < ^ (21B) 

regressive model. However, in order to validate the recog- 2jN 2MinPower - 

nition and arrive at a conclusive detection, preferably some s 

additional decisional criteria are added. This added process- , 

inn is preferred due to the presence of noise and distortion where Mmpower is the threshold value 

in the analyzed signal as a result of which the a 2 coefficieut u Conversely, any of the fo lowing conditions may be met 

will not exactly match the value of p In addition, the value before decknn S *° termination of a single frequency tone: 

of the coefficient a, depends only on the decay factor p and 10 a ) not «>°«g»» P° wer 1Q ,he aaalyzed signal, e.g., 

the frequency of the tone. The analyzed signal is considered <Mlnr °" er > . 

to be sinusoidal only when the modeled signal has the •>),*» •? coefficient is not close to p , i.e., (p 

correct decay factor. In this case, the coefficient a a should <P an f optonafly, 

depend only on the frequency (since the decay factor is <=) the a. coefficient is notstable, i.e., to «e one or more 

known and fixed) but due to the presence of noise or 15 a ' ( n > wth <*i(0-°W«») or "■(■Wa 1 (#) within 

distortion,therealvalueofa 1 willbeinfluencedbythevalue Tte'firsfcondition preferably causes the detector 10 to 

°* a 2 > , _ , declare the loss of the tone at any sample where the estimate 

Therefore, to validate the presence of a tone, some degree of me analyzed signal power becomes less than the 

of debouncing may be required. For instance, the value of minimum required (Minpower). The second condition 

the auto-regressive coefficient 2^ should preferably be equal 20 snou ^ oe considered true when the condition is met over a 

or close to p 2 within an acceptable tolerance e over a number of N off consecutive samples. By requiring consecu- 

contiguous number of samples before the signal is declared ti ve samples to meet the condition, the detector is made 

to be sinusoidal in order to take in account noise and more immune to noise and distortion, 

distortion. The condition of the a^ coefficient being close to Note that the decay factor p preferably lies in the range of 

P 2 can be mathematically expressed as: 25 0<p<l. When p is a relatively low value, the exponentially 

decaying sine-wave equivalent of the signal under analysis 

(p 2 -€)<a 2 <(p 2 +t) (19) retains a weak "memory" of prior samples of the signals. 

Conversely, when p is a relatively high value, the exponen- 

and preferably this condition will exist over Inconsecutive ^ decaying sme . waV e equivalent of the signal under 

samples before declaring the presence of a single frequency 3Q mdlysis retains a strong "memory" of prior samples of the 

tone. Alternatively, the value of the coefficient^ may be sign ^ j^s characteristic can be employed to tune the 

within some other predetermined range of p over N 0 „ detector. More specifically, when P is a relatively high value, 

consecutive samples. tne speed of detection is decreased but the accuracy is 

If desired, an optional extra condition to use in the i ncrea sed, and when P is a relatively low value the speed of 

validation is to test for an almost constant value for the 35 detection is increased but the accuracy is decreased, 

coefficient a lf thereby indicating a tone with stable fre- FIG. IB is a flow chart summarizing the above discussion 

quency. This condition can be mathematically expressed as: and demonstrating the program 20 executed on the DSP 12. 

, , A ^ fcr • , MA In a first step 22 all the required variables are initialized. 

MO-^^MO**), for «^«<^i (20) specifi £ aUy> me aut0 _ regressive coefficients a, and a, 

where a^n) is the current a, coefficient, a a (i) is one of a 40 are set to zero, the correlation coefficients a', 0 ,a' a' 0 ,a' 12 , 
number of previously calculated a, coefficients, and 8 is a <x' 21 ,a' ?2 are set to zero, and a three-deep sample history 
tolerance error. In the preferred embodiment, the value of a, b ^ 10 zeio v 

is computed from a, in accordance with equation (17A) or Step 24 through 38 form part of an endless loop construct. 
(17B) and averaged over N w samples in order to report the At step 24, the signal 14 is sampled at a rate f, and the 
value of the detected frequency. 45 sample history buffer is updated. 

Another preferred criteria in validating the presence of a At a following step 26, the correlation coefficients a 10 , 
tone is to ensure that the power of the incoming signal is fnA^A^ are computed .preferably in accor- 
sufficiently high to ensure non-spurious readings. The coef- dance with one of the equations sets (10) or (18). Note that 
ficient a 22 is readily available from the above correlatioo the correlation coefficients are updated for each sample of 
coefficients and provides a good approximation of the 50 ^ sl S nal obtained at step 24. 

incoming power if the decay factor a is close enough to At a following step 28 a counter is tested. The counter 
unity. This is because under such circumstances a 22 is determines whether or not P samples have been obtained. If 
orooortional to contro1 B passed to the folio wing step 30 for compu- 

tation of the auto-regressive coefficients and subsequent 
55 decision of whether or not a tone has been recognized; 

ZS t otherwise control is passed back to sampling step 24. The 

N value of P may be 1, in which case steps 30 and following 

are performed for each sample of the signal. Alternatively, 
the power of the signal, where N is the number of samples these steps may be periodically performed in which case P 
used to perform the power measurement. Alternatively, any 60 will be greater than 1. The choice of the value of P will 
other measure of the power of the incoming signal can be depend on the specific application and the sampling fre- 
used, including a running computation of the foregoing quency. 

equation. Thus, one way the condition of the existence of At step 30 the counter used in step 28 is reset and the 
sufficient power in the analyzed signal can be mathemati- auto-regressive coefficients a A and a^ are computed in accor- 
cally expressed as: " 65 dance with equations (16) and (17). 

At step 32 the tone recognition condition specified in 
022^ MinPowcr, (21A) equation (19) is evaluated. Preferably the conditions speci- 
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fied in equations (20) and (21) are also evaluated in order to of the square of an exponential decay factor used to define 

validate the tone recognition, and the non-fulfilment of any the exponentially weighted future sliding equivalent of the 

of these conditions may indicate the absence of a tone where signal. 

presently one had been detected. 4. The method according to claim 3, wherein the step of 

At step 34, if the conditions evaluated at step 32 have been s recognizing the presence of the tone further includes vali- 

fulfilled then control passes to step 36 which sets the dating the presence of the tone by determining the power of 

detection output signal 16 and optionally the numerical the signal and confirming that the power exceeds a threshold 

value of the detected frequency as evaluated by equation power. 

(15) is provided on frequency output signal 18. Otherwise, 5 memo d according to claim 4, wherein validating 
control passes to step 38 wherein the detection output signal 10 ^ preserice 0 f me tone further includes computing the first 
16 is reset if previously set or left unchanged. auto-regressive coefficient and confirming the stability of the 
FIG. 2 shows how the tone detectors 10 may be used in first au to -regressive coefficient over a number of iterations, 
an echo-cancellation system 40. Two tone detectors 10 are 6 ^ met hod according to claim 1, wherein the auto- 
disposed to intercept bi-directional signals propagated regressive model for the exponentially weighted fixture 
through an adaptive echo-cancellation circuit 42, as shown. 15 sliding equivalent of the signal is defined by a transfer 
The detection and frequency outputs 16 and 18 of the tone function H(z), wherein: 
detectors are connected to a decision logic circuit 44, which 

in rum is connected to the echo-cancellation circuit 42 so as A 

to freeze or disable it. The logic circuit 44 is configured to = (2xfo) ~ 

disable the echo-cancellation circuit 42 if any of the tone 20 1 -2/fcosJ-^-Jr 1 i-pr 1 
detectors detect a tone of a specific frequency used to convey 
signalling information in order to avoid the problem of 

eliminating the tone as described above. The logic circuit 44 where fo is the frequency of a detected tone, 

is also configured to freeze, i.e., stop the adaptive behaviour f, is the sampling frequency, and 

of, the echo-cancellation circuit 42 if both tone detectors ^ p is an exponential decay factor. 

detect a tone of the same frequency. 7. The method according to claim 6 wherein the first and 

Those skilled in the art will appreciate that numerous second auto-regressive coefficients, a 3 and a^ are defined 

modifications and variations may be made to the preferred by: 
embodiment without departing from the spirit and scope of 

the invention. 30 

What is claimed is: 

1. A method for recognizing a tone in a signal using at 
least a second order non-biased auto -regressive model a 2 =J?- 
defined by a first auto-regressive coefficient and a second 

auto -regressive coefficient, said method comprising: 35 

(a) sampling the signal at a constant sampling rate and, for 8 The metnod according to claim 7, wherein the current 
each sample thereof, recursively determining a finite correlation coefficients are d 10> a\ Xt a.\ Oy a\^ lt a: 22 
number of current correlation coefficients usmg an wfaich m computed according to one of the following first 
exponentially weighted, future sliding equivalent of the ^ second ^ of equations> the first set being> 

signal, said signal being time reversed such that prior to 40 

each determination of the current correlation coeffi- a'^p'x^+p 2 ^,, 

cients as aforesaid, a current received sample and 22 

consecutively received previous samples of the signal a ' 20 "^ ^ ^ 

are time reversed such that the said current received a'^-f^x^ 2 ^^ 

sample is defined as an oldest sample, an oldest of the 45 

said consecutively received previous samples is defined a'i2-r^i*o+P 2 "i2 

as a current sample and, where more than two con- qI + ^2 a 

secutively received previous samples of the signal are 21 1 21 

used in the said determination, all said consecutively a^-jco^+P 2 ^ 

received previous samples other than the said oldest 50 

thereof are also denned to be time reversed as to their tnc second sct bein S> 

received sampling order, and wherein the current cor- a' 10 =p 2 a 12 

relation coefficients are determined using pre-existing 

values of correlation coefficients determined in a pre- a'n=P 2a 22 

vious determination, the current sample and at least two 55 a^^jto^+p 2 ^ 

consecutively received previous samples of the signal; 

(b) periodically determining at least a second auto- a'i2=fco*i + P 2a i2 
regressive coefficient for modelling the signal using the a^^+^a^ 
current correlation coefficients; and 

(c) recognizing the presence of the tone based on the 60 a 2i- a 'i2 
second auto-regressive coefficient. wherein* 

2. The method according to claim 1, wherein step 1(6) is 

executed every P samples, P being an integer number of at *o * the current sample of the signal, 

least one. x i is a sample of the signal immediately previous to x 0 , 

3. The method according to claim 2, wherein recognizing 65 Xj is a sample of the signal immediately previous to x Jt 
the presence of the tone includes confirming that the second and a 10 ,a 11 ,a 20 ,a 12 ,a 21 ,a 22 are the pre-existing values 
auto -regressive coefficient is within a predetermined range of the correlation coefficients. 
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9. The method according to claim 8, wherein the first and sliding equivalent of the signal is defined by a transfer 
second auto-regressive coefficients a x and a 2 are computed function H(z), wherein 
from the current correlation coefficients as follows: 



10. A method for detecting a tone and its frequency in a 
signal using at least a second order non-biased auto- 
regressive model defined by a first auto-regressive coeffi- 
cient and a second auto-regressive coefficient, said method 
comprising: 

(a) sampling the signal at a constant sampling rate and, for 
each sample thereof, recursively computing a finite 
number of current correlation coefficients using a an 
exponentially weighted future sliding equivalent of the 
signal, said signal being time reversed such that prior to 
each computation of the current correlation coefficients 
as aforesaid, a current received sample and consecu- 
tively received previous samples of the signal are time 
reversed such that the said current received sample is 
defined as an oldest sample, an oldest of said consecu- 
tively received previous samples is defined as a current 
sample and, where more than two consecutively 
received previous samples of the signals are used in the 
said computation, all said consecutively received pre- 
vious samples other than the said oldest thereof are also 
defined to be reversed as to their received sampling 
order, and wherein the current correlation coefficients 
are computed using correlation coefficients computed 
in a previous iteration, the current sample of the signal 
and at least two consecutively received previous 
samples of the signal; 

(b) periodically computing at least the first auto- 
regressive coefficient and the second auto-regressive 
coefficient for modelling the signal using the current 
correlation coefficients; 

(c) recognizing the presence of the tone based on the 
second auto -regressive coefficient; and 

(d) determining the frequency of the tone based on the 
first auto-regressive coefficient and the sampling fre- 
quency. 

11. The method according to claim 10, wherein step 10(&) 
is executed every P samples, P being an integer number of 
at least one. 

12. The method according to claim 11, wherein the step of 
recognizing the presence of the tone includes confirming 
that the second auto -regressive coefficient is within a pre- 
determined range of the square of an exponential decay 
factor used to define the exponentially weighted future 
sliding equivalent of the signal. 

13. The method according to claim 12, wherein the step 
of recognizing the presence of the tone further includes 
validating the presence of the tone by determining the power 
of the signal and confirming that the power exceeds a 
threshold power. 

14. The method according to claim 13, wherein validating 
the presence of the tone further includes confirming the 
stability of the first auto-regressive coefficient over a number 
of iterations. 

15. The method according to claim 10, wherein the 
auto -regressive model for the exponentially weighted future 



where is the frequency of a detected tone, 
10 f, is the sampling frequency, and 
p is an exponential decay factor. 

16. The method according to claim 15, wherein the first 
and second auto-regressive coefficients, 2l x and ^ are 
defined by: 

15 

ai = and 

00 

and 

17. The method according to claim 16, wherein the 
current correlation coefficients are &\ Q ,a\ li o! 2Q9 a? vl> a! 71t 
a! 12 which are computed according to one of the following 

25 first and second sets of equations, the first set being, 

and the second set being, 

Ctll=P <*22 

a'2o=P 2 ^2 + P 2 ci2o 

<*'l2-fo>*l+P 2 <*12 

a, 2t B,a 'i2 
wherein: 

50 Xq is the current sample of the signal, 

x-l is a sample of the signal immediately previous to x 0 , 
^ is a sample of the signal immediately previous to x^ 
and a 10 ,ct 1:l ,0120,042,012! ,a 22 are values of the correla- 
tion coefficients computed in a previous iteration. 
55 18. The method according to claim 17, wherein the first 
and second auto -regressive coefficients a x and a2 are com- 
puted from the current correlation coefficients as follows: 

<*21<*io 1<*20 , „ tt 10 QT 12 a 2 

An a-2 = , and a y = . 

***** 
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